Transport anomaly in the low energy regime of spin chains 
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The anomalous thermal conductivity in spin chains observed in experiments is studied for the 
low temperature regime. In the effective dynamics with most realistic perturbations, the so-called 
Umklapp terms is irrelevant to reduce mean free path in the energy transport at even finite temper- 
atures. This is consistent with large conductivities found in recent experiments. The Drude weight 
which is the prefactor in the divergent conductivity is calculated, and the temperature dependence 
is discussed. 
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I. INTRODUCTION 

Recent experiments on transport phenomena in 
low dim^nAkmal^ipin systems have attracted much 
interestEI^'LKffl'LHjij. There the behavior of conductiv- 
ity is sensitively affected by the magnetic phases. There- 
fore the thermal conductivity can be a strong experimen- 
tal instrument to investigate the magnetic properties in 
materials. In experiments, magnetic energy transport is 
much larger than the phononic transport at most tem- 
peratures. This fact is considered to be attributed to 
large exchange coupling constants and nondiffusive be- 
havior of transport. Sr 2 Cu03 described by the isotropic 
Heisenberg chain has large exchange coupling constant 
J 2000-R' and the magnon can transmit energy even at 
the order of 100K. The mean free path is experimentally 
estimated in the similar way as JJiCL-Pcicrls-Boltzmann 
approach for phononic transportli!JO, and nondiffusive 
behavior is found. In another isotropic Heisenberg chain 
CuGeC>3, where the mean free path is calculated as 500 
times lattice size. It is noteworthy that even when the 
alternation between bonds occurs, the magnetic trans- 
port may be still dominanLand nontrivial enhancement of 
conductivity is observedtrtj. Spin-ladder system (Sr, La, 
Ca)i4Cu2404i and L|%i < fff*pC'U2404i also show extremely 
large mean free pathQU'EfLj. 

As studied in classical systems^ an integrable system 
shows ballistic energy transportlij. Even if the energy 
current is not conserved, the low-dimensional systems 
with conserved quantities tend to show the divergence of 
conductivityEj. . Xhrae rtepdeBfiies are also valid in quan- 
tum systemMByOBBBy. The isotropic Heisen- 
berg chain is the integrable system and exactly shaws the 
diverging conductivity as proved by Zotos et al.EJ. This 
divergence is a direct result of conserved energy current. 
However in most materials in experiments, the perturba- 
tions like a bond-alternation and a next nearest neigh- 
bor interaction exist, so that the energy current is no 
more conserved. Nevertheless the experiments indicate 
divergent conductivities. Motivated by these experimen- 
tal facts, Alvarez and Gros investigated whether the ex- 
perimentally indicated divergent conductivity really oc- 
curs or not. They numerically found nonvanishing Drude 



weight at finite low temperatures for the alternate spin 
chain with frustration and the ladder systemEEl, although 
it may vanish, at very high temperatures in the thermody- 
namic limittj. Numerical finding of finite Drude weight 
indicates the existence of divergent conductivity at low 
temperaure regime. Stimulated by these numerical stud- 
ies, in this paper we study the possible mechanism to 
explain nondiffusive behavior in the presence of perturba- 
tions to-the Heisenberg chain. We focus on the low energy 



regimecrl by employing the effective dynamics as adopted 
for the Heisenberg chain by Meisner et al.Ea. As a result, 
we find the effective energy current becomes conserved 
in spite of the existence of so-called Umklapp terms so 
that the conductivity diverges. This ballistic nature is 
a common characteristic for most realistic perturbations 
as far as the low energy regime is considered. This is a 
remarkable aspect in the sense that the perturbations is 
irrelevant to decrease the mean free path although they 
can be relevant to cause various magnetic phases in the 
equilibrium state. This is very important to explain the 
anomalous large conductivities of one-dimensional mag- 
netic systems. We further calculated the quantitative 
value of the Drude weight of conductivity using the path 
integral approach, and found it to be finite. 



II. EFFECTIVE ENERGY CURRENT IN LOW 
TEMPERATURE REGIME 

We first consider an alternate spin chain whose Hamil- 
tonian is 

H = H +H U (1) 

H = jJ2SfS*+i, Hi = J<5^(-l) £+ %-S £+1 . (2) 

t i 

The total Hamiltonian is divided into the isotropic part 
7io and the dimerized part Hi- Thermal conductivity 
k{T) at a temperature T(= 1//3) is, calculated by the 
Green-Kubo formula which reads asE3, 



k(T) = Urn ik(l>) = -j^- 



die; 



dT(j(t-ir)j), (3) 
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where is the total current operator in the Heisenberg 
picture. From the continuity equation of energy, the total 
current operator j is written as j = —i J2e ht+{\ with 
the local Hamiltonian hi := J(l + {-Vf +1 8)Si ■ Sf +1 Ei 
In the absence of the alternation, j is the conserved 
quantityEJ, whereas the existence of alternation does not 
allow the conservation, i.e., 



[i(a=o).Wo] =0, but \j,H]^0. 



(4) 



At low energy regime, the Hamiltooian can be effec- 
tively represented by the boson fieldc^tJ; 



dx 



irvp 2 + —(d x e+(x)) 2 - BcosO + (x) , (5) 



4/T 



where v — 4? and B = J5. The Boson field is defined 
using the fermion field and and satisfies the com- 
mutation relation as, 



I dB± 
2ir dx 

p(x) 



tf R {x)9 R (x):±:9l(x)9 L (x) :, (6) 
1 



4tt 



d x Q-(x) 



[0 + (x),e_(x')] = 27risgn(x- x'). 
We define the local Hamiltonian density at x as, h(x) ■ 
nvp 2 + \ ~ £?cos0+(x). Using the relations, 

56+ 



(7) 

(8) 



-2nvp, 



dp 



V 



d x 2& + {x) + 5sin0+, (9) 



dt " Jrl dt 2tt' 

we find the Heisenberg equation of motion given by, 



dh(x,t) v* 



d x ( P d x Q+ + (d x e + )p). 



(io) 



dt 2 

This leads the effective total energy current written as, 



3b =v 2 dxp(x)d x Q + (x). 



(11) 



Interestingly, this effective current operator j'b docs 
not include the perturbation term (Umklapp term) 
-Bcos0+(a;). Furthermore we find that even if the en- 
ergy current is not rigorously conserved (^), it behaves 
as the conserved quantity in the low temperature regime, 
i.e., 



\jH] ^0, but [j B , H B }=0. 



(12) 



The perturbation term £?cos0+(a;) causeSpthe Umklapp 
process in the scattering between fermion£3. Neverthe- 
less at the level of the energy transport, this process is 
irrelevant to decrease the mean free path at finite low 
temperatures. 

This situation is also the case in the other realistiCrper- 
turbations. We consider the two following examplesEa, 

n 2 = Ja^s,-s, +2 , n 3 = j 7 ^(-r/sf, (13) 



The perturbation 7i 2 causes the frustration, and the non- 
trivial phases appears as a function of o£3. Recent nu- 
merical studies treat this .pjEturbation and the effect of 
frustration is investigatecESlia. H3 is the term of stag- 
gered magnetic field. The effective expressions at the low 
energy regime for the perturbations ([L3|) are written as, 

H 2 ~ J dx c 2 (d x <d+f + c' 2 cos (29+) , (14) 
T-i 3 ~ c 3 J dx cos (26+) , (15) 

with constants, c 2 , c' 2 and C3. These cases have the same 
form as the Hamiltonian (H), and the effective total cur- 
rent is easily proved to be conserved. This feature is 
generally common for the Hamiltonian with the form, 



dx 



ap 2 + b{d x Q+{x)) 2 



ci cos (diQ+{x) 



(16) 



Here a,b,ce,di, and (j>i are constants. Thus we conclude 
that the effective energy current is conserved in most re- 
alistic spin chains at least as far as we consider in the low 
temperature regime. 



III. DRUDE WEIGHT 

Now let us calculate the prefactor of thermal conduc- 
tivity, i.e., the Drude weight Dth in the low energy regime 
and confirm that finite value of Dth really exists. Since 
the energy current does not depend on a time, we write 
the Green-Kubo formula as, 

k = D th S{u) (17) 

D th =/3Tr/L[ dr f dx I dx' (j B (x,T)jB{x',0)), (18) 
Jo Jo Jo 

where j(x,r) is the local energy current at x and the 
imaginary time t in the Heisenberg picture. The current 
is expressed as j(x,r) = d x <d + (x, t) d T Q+(x, t). We 
confine ourselves in the alternate spin chain effectively 
described by Hamiltonian (^) . At very low temperatures, 
the Haaiiltonian (|5|) is approximated by the mean field 
theoryE3 ; 



cos(e+)~e-< e +>/ 2 l-^± 



e 2 + - (el) 



(19) 



The Hamiltonian is reduced to 
Hb ~ / dx 



— (d x e+(x)) 2 + C9i + nvp 2 + const. 

47T 



(20) 



where C — Be ( e +^ 2 and (0+) is the average value at 
the ground state which is self-consistently determined. 
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This approximation yields the accurate ground state en- 
ergy within 10 percept error in comparison with the ex- 
act numerical datac3. Thus we expect that this approx- 
imation well describes properties in very low temper- 
atures. Using the Hamiltonian (p0|), we calculate the 
Drude weight D t h using the path integral (e.g.,E£l). The 
Drude weight is given by explicit calculation of 



D 



tb 



4tt/3L 3 



J J J drdxdx' d Xl d X3 d T2 d Ti ^ 



ni,..,?i4 ki,..,k4 



x e < i{kiXi J r k2X 2 -\-k 3 x 3 -\-\-k J iX4) ^i(uj ril ri+^,i 2 r 2 +aj Tl3 r 3 ++aj Tl4 T 4 ) 

x (e + (*i,ni)e + (fe,na)e+(fc3,n3)e + (fc4,n4)), (21) 

where we take x\ = X2 = x, X3 — £4 = x' , t\ = T2 = r, 
andr3 = T4 = 0. k m and uj n are the wave number 
and the Matsubara frequency defined as k m — and 
u n = ^jp, respectively. 9+(fc m ,n) is the Fourier trans- 
formation at the mode k m and w n . The four point cor- 
relation is calculated by the Wick's theorem using the 
two point correlation (|8+(fc m , n)\ 2 ) = ^ +v ^ +i7rvC ■ 
We next carry out the exact summation over n p (p — 
1, • • • , 4) using technique of the complex integral. We fi- 
nally obtain the simple expression of the Drude weight 
with a — Jr-\/v 2 k 2 + 4nvC, 



1 e- 



1 e- 



2vT / dx , 

lo sinlr yjx 1 + ir/3 2 vC 



(22) 



In the case of isotropic Heisenberg chain C = 0, Dth = 
VV 3 T because of J °° dx ^ = \ , which completely 
agrees with the asyntatotic behavior of exact calculation 
by the Bethe ansatzcil. This .result is also given by the 
conformal transformationEJ'Eil. In the presence of finite 
alternation S, however, the Hamiltonian does not have 
conformal symmetry 

Eq. ( p2| ) shows that the Drude weight is non- vanishing 
for finite alternation S. This is quite crucial because the 
existence of Dth assures the anomalous energy trans- 
port. The energy gap AE is roughly estimated as 
AE ~ 2\ttvC. Thus for very low temperature regime 
T <C AE, eq.(p2|) shows exponential decrease with power 
law prefactor due to the finite gap 



D 



th 



Vf3 2 f°° , 2 -Ps+Aj 



v 
32 



ttAD 3 



-0 AE 



T 



T < AE. 



(23) 



Exponential decrease is simply due to exponential de- 
crease of thermal excitation. The power law prefactor 



T- 1 ' 2 is also seen in the electric aad-apin conductivity 
in one-dimensional gapped systemsf3Ej. This tempera- 
ture dependence will be the characteristic in transport 
phenomena in the low-dimensional gapped systems. 



IV. SUMMARY 

We find the possible mechanism of anomalous trans- 
port in low energy region of realistic spin chains. In the 
isotropic Heisenberg chain, the energy current is exactly 
conserved so that the Green-Kubo formula trivially di- 
verges. On the other hand, the presence of perturbations, 
the Umklapp terms appear. However this term is irrel- 
evant to the energy current of long wave length at low 
temperatures. We considered the quantitative estimation 
of the Drude weight for the alternate chain, and find ex- 
ponential decrease of Drude weight at low temperature. 
This aspect is consistent with the numerical studiesE3E3. 

Here we should make clear the difference between claa=. 
sical nonintegrable Fermi-Pasta-Ulam (FPU) chainsllj 
where the mechanism of anomalous transport is still ar- 
gued. In this classical system, anomalous heat transport 
is attributed to a slow relaxation due to the total momen- 
tum conservation, whereas the the auto-correlation func- 
tion of energy current vanishes in the long time limit, 
i.e., the mixing property is satisfied. However, in per- 
turbed quantum Heisenberg spin chains we disewfis here, 
Numerical evidence of the finite Drude weighiil3 means 
the violation of the mixing property. Thus the mecha- 
nism in quantum spin chains is different from the one in 
classical FPU system. 

The alternate spin chain directly related to the spin- 
Peierls system. In the realistic mechanism, the alterna- 
tion is formed by the fluctuation of the phonon. In this 
sense, in more realistic treatment, 5 must be exchanged 
by the time-dependent bond length between atoms. In 
the case where the bond fluctuation is not negligible, how 
the magnetic energy transport is affected is also interest- 
ing. 

note added- After submitting this paper, we find the 
preprint! 3 -! where the Drude weight in the dimerized xy 
spin chain and two-leg ladder system are extensively 
studied. The paper includes the similar formula as (122 
in these systems. 
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